Two-mode Bose gas: Beyond classical squeezing 
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The dynamical evolution of squeezing correlations in an ultracold Bose-Einstein distributed across 
two modes is investigated theoretically in the framework of the Bose-Hubbard model. It is shown 
that the eigenstates of the Hamiltonian do not exploit the full region allowed by Heisenberg's un- 
certainty relation for number and phase fluctuations. The development of non-classical correlations 
and relative number squeezing is studied in the transition from the Josephson to the Fock regime. 
Comparing the full quantum evolution with classical statistical simulations allows to identify quan- 
tum aspects of the squeezing formation. In the quantum regime, the measurement of squeezing 
allows to distinguish even and odd total particle numbers. 

PACS numbers: 03.67.Bg, 03.75.-b, 03.75.Gg, 03.75.Lm, 05.20.-y, 05.30.-d HD-THEP-10-04 



I. INTRODUCTION 

Precision measurements at the quantum level are ulti- 
mately limited by Heisenberg's uncertainty relation. De- 
creasing the fluctuations of an observable of interest be- 
low the standard quantum limit given by the central limit 
theorem necessarily increases the fluctuations in one or 
more conjugate observables. Different ways to use such 
squeezed states, e.g., to measure frequency in Ramsey- 
type interferometers have been discussed in great detail 
in the seminal papers [ll-Q- Heisenberg- limited Mach- 
Zchndcr intcrfcrometry using number squeezed photon 
states were studied in detail in . The standard quan- 
tum limit is reached by today's best sensors of various 
quantities such as time @ and position (Tol [Tl| . 

Two-mode Bosc-Einstcin condensates of non inter- 
acting particles constitute coherent semiclassical macro- 
scopic ensembles of particles. In the limit of zero tem- 
perature the probability for each mode to contain a cer- 
tain number of atoms is approximately Poissonian. The 
variance (An) 2 of the particle number difference is pro- 
portional to the total particle number N, corresponding 
to the standard quantum limit (An)/N ~ \/y/N. Con- 
jugate to particle number difference is the relative phase 
4> between the modes which can be measured through 
interference effects [49[. The Heisenberg uncertainty re- 
lation AnA</> ~ 1 implies the relative standard deviation 

(A0)/0 ~ i/Vn. 

The presence of interactions between the particles 
strongly modifies the situation. Fluctuations of the par- 
ticle number and squeezing in trapped atomic gases have 
been the subject of numerous recent experimental and 
theoretical studies Here we are interested in 

squeezed states of Bose-Einstein condensed ensembles of 
trapped atoms. We specifically consider squeezing in the 
particle number difference between the two minima of a 
double-well trap in one spatial dimension. This reduc- 
tion of number fluctuations below the standard quantum 







Sque< 
Me 

■ 111 


azing 

trology gain 



0.2 0.4 0.6 0.8 1 



\a\ 

FIG. 1: (Color online) Phase diagram for N = 100 particles 
in a double-well potential, with on the average equal popula- 
tions in the two wells, (h) = {hi — U2)/2 = 0, and the abso- 
lute phase chosen such that {§2) = i{a\a\ — a\a,2) = 0, see 
Eq. (|A9[) . For any given coherence parameter a = 2{Si)/N 
the variance (An) 2 = of the number difference be- 

tween the wells is bounded below and above by Heisenberg's 
uncertainty relation and the constraint that the total parti- 
cle number N is fixed. The allowed area is defined by the 
inequality (|11[) shaded in colors in the above diagram. In the 
medium (violet) and dark (blue) shaded regimes, the squeez- 
ing parameter is below the standard quantum limit, £3 < 1, 
see Eq. ([5]). In the dark (blue) regime, the squeezing can be 
used to gain precision in metrology, see Eq. (TH)) . 



limit occurs at the expense of increased fluctuations in 
the relative phase between the wells. The barrier be- 
tween the wells of the system is taken much smaller than 
the outer walls of the trap such that the main source of 
number fluctuations in each well is given by tunneling 
processes between the wells. Such two-mode quantum 
systems can be described, by use of Schwinger's repre- 
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sentation, in terms of angular momentum states with the 
maximum length of the spin vector related to the total 
particle number. Consequently, the non-classical states 
we consider exhibit a variant of spin squeezing P, Q ■ 

Our studies have been initiated by a recent experiment 
[ijj in which a setup of the kind sketched above has been 
used, with a one-dimensional double- well trap created by 
a superposition of standing light waves. In this experi- 
ment, both particle number and phase difference have 
been measured independently with high resolution to es- 
tablish squeezing correlations empirically. While our dis- 
cussion of the formation of squeezed states here focuses 
on particle number variations in space, it can be straight- 
forwardly extended to apply to alternative schemes, in- 
cluding, e.g., squeezing in the relative occupation number 
of internal hyperfine states of the atoms [19( . 

A prominent motivation of the work presented here is 
to study the role of quantum statistical fluctuations in 
the preparation of the spin squeezed states. By compar- 
ison of exact quantum with semi-classical Monte Carlo 
simulations it is shown here that the spin s quee zing pro- 
duced in a setup as that in the experiment [181 ] is equiv- 
alent to reduced classical fluctuation of the occupation 
number difference between the two wells. It is, further- 
more shown that the production of squeezed states in the 
experiment follows a quasi static path in state space. 

Spin s quee zing is closely related to quantum entangle- 
ment [2CH26|. Schemes have been proposed to use spin- 
squeezed states for quantum t elep ortation of continuous- 
variable, e.g. coherent states |27H29| . However, a full re- 
construction of the teleported state is only possible in the 
limit of perfect squeezing (30j . close to zero temperature 
where quantum fluctuations become relevant, beyond the 
regime of classical squeezing. We demonstrate to which 
extent the evolution of the system enters the regime of 
non-classical states. 

We furthermore show that, for a finite total number 
of atoms the Heisenberg limit can not be reached, in the 
generic case, by adiabatically changing the system's pa- 
rameters, starting from an incoherent mixture of energy 
eigenstates. In particular, in the ground state of the sys- 
tem, though being rather close to the Heisenberg limit, 
the variances of the conjugate variables number differ- 
ence and relative phase between the wells are in general 
not minimized at the Heisenberg limit. The difference 
between ground state and Heisenberg limited variances 
becomes negligible only in the limit of large total parti- 
cle number. Any excited energy eigenstate is found to 
be even further from this limit. The production of states 
at the Heisenberg limit therefore requires nonequilibrium 
dynamical evolution. 

Our paper is organized as follows: In Sect.|H]we discuss 
fundamental squeezing limits in the Bloch-sphere picture. 
The dynamics of the production of squeezing by raising 
the barrier in the double-well potential is described in 
Sect. IIIII A focus is set on comparing the full quantum 
with the semiclassical evolutions. Our conclusions are 
drawn in Sect. IIV1 



II. NUMBER-SQUEEZING IN A TWO-MODE 
BOSE GAS 

We consider a Bose condensate trapped in a double- 
well potential with two energetically degenerate minima 
separated by a barrier of variable height, allowing for 
an adjustable tunneling rate between the wells. In the 
experiment (l8j such a potential was formed optically by 
counterpropagating laser waves. We consider parameter 
regimes in which the system can be described by the two- 
site Bose-Hubbard Hamiltonian 

U 2 

H = - J(a\a 2 + d|oi) + Tj- ^2 ajajoifij. (1) 

i=l 

where the operators a, and a] obey the standard bosonic 
commutation relations 

[ai,a,j]=0, [ot,ot]=0, [ai,aj} = Sij. (2) 

and where J is the tunneling and U is the onsite interac- 
tion parameter. The energy spectrum and eigenstates of 
this Hamiltonian, as far as relevant for the discussion in 
this article, are summarized in Appendix [A] The Hamil- 
tonian can be written in terms of Schwinger angular mo- 
mentum operators, 

H = —J(S + + S-) + U (Sf + S 2 — N), (3) 

where S± = Si ± iS 2 , and S 2 = J^Li with Si = 
{a\a\ + ajo,2)/2, §2 — i{a\ai — d{a2)/2, and S3 — 
(a[ai — a 2 a2)/2. The corresponding angular-momentum- 
type states of the system can be represented by quantum 
phase-space (Wigner) distributions on the Bloch sphere 
[471 [48| . For details about this representation we defer 
to Appendices [5] and [B] 

The difference in occupation number between the 
modes is measured by S3, while the relative phase 
coherence is measured by the orthogonal components 
Si^2- These angular-momcntum-type obscrvables are 
constraint by a Heisenberg uncertainty relation, see 
Eq. (|8|). which allows us to derive fundamental limits for 
squeezing. 

To quantify angular momentum squeezing one intro- 
duces a squeezing parameter adapted to the problem 
under consideration. A definition suitable for number 
squeezing between bosonic modes relates the variance of 
one Cartesian component to the total spin S = N/2, i.e., 
to the total number of mode excitations or particles, 

& = (ASrf/iS/2). (4) 

For z = 3, this is equivalent to the ratio of the variance of 
the number difference m — n 2 — 2n to the total number, 

£3 = 4(An) 2 /A. (5) 

In the following section, we show the evolution of this 
squeezing parameter under the slow ramp-up of the po- 
tential barrier. It essentially reflects the suppression of 
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FIG. 2: (Color online) Positions of the energy eigenstates 
of the Hamiltonian Q for a system with N = 100 particles 
in the phase diagram introduced in Fig. [T] Each black dot 
corresponds to one such state. The dotted lines are drawn to 
guide the eye between successive eigenstates of a Hamiltonian 
with fixed ratio U/J. From top to bottom, the sets of states 
are obtained for U/J = 0,1,10, and 100, respectively. The 
(red) solid line connects all ground states for different U /J. 
We emphasise that the energy eigenstates do not extend over 
the full (shaded) region allowed by Heisenberg's uncertainty. 



the relative number fluctuations below the classical limit 
(An) 2 = N/4 as given by the central-limit theorem. 

In the following we will consider only situations where 
the total number TV of atoms in the wells and therefore 
the total spin S are fixed. Before proceeding with this we 
briefly remark that in the general case of varying S the 
squeezing parameters £j measure the deviation from an 
angular momentum coherent state which is represented, 
in Bloch space, as a spherical Gaussian uncertainty dis- 
tribution around (S) with radial width a = ^f~N /2. Such 
a state can be written as a product of coherent states in 
the two modes 1 and 2, |a)|a), |a| = \J iV/2, and has 

ii = 6 = £ 3 = 1. 

In analogy to spin-squeezed states, two-mode states 
can exhibit squeezing in particular directions at the ex- 
pense of increased fluctuations in directions perpendicu- 
lar to this 0]. For example, an eigenstate of S3 can be 
illustrated by a circle on the Bloch sphere parallel to the 
1-2-planc corresponding to the variances 



(ASr) 2 



(AS 2 ) 2 
(AS3) 2 



= [S(S 

= 0, 



1) 



(6) 



respectively The operators Si , S2 measure the phase be- 
tween the Fock modes 1 and 2 which is accessible to in- 
terference measurements of the particle occupation num- 
bers, see Ref. Hence, the commutators (|A10[) char- 
acterise the uncertainty relation ([5]) between (relative) 
number and cosine of the (relative) phase. The absolute 
phase is undefined as we assume a fixed total particle 



number. Physically this phase is not measurable without 
comparing and therefore coupling the system to another 
system. Hence, the two-mode states considered in the 
following, with relative particle number centered around 
n = 0, correspond to a distribution on the surface of the 
Bloch sphere, centered around the equator, i.e., around 
a polar angle 8 = it/2 or (S3) = 0. A more convenient 
parameter to reflect the squeezing below the standard 
quantum limit, in accordance with the uncertainty rela- 
tion (|SJ) is given by 



£« - (AS fc ) 2 /|(5)/2|. 



(7) 



In the following, we will consider the special case (S2) = 
(S3) = such that the only nontrivial combinations are 
£21 and £31. The variances (ASfc) 2 are subject to the 
uncertainty relation 



(AS fc ) 2 (AS z ) 2 > \\e Um {S m )\ 2 - 



(8) 



such that the parameters £&; obey the inequality 

£ 2 i£ 3 i > 1- (9) 

This defines the Heisenberg limit for the fluctuations. 
One may at first glance expect that the limit £31 — > 
is allowed at the expense of £21 — ► °o and vice versa. 
However, as long as |(Si)| > perfect squeezing in ei- 
ther the relative number or the phase is not possible as 
we have assumed fixed S and thus a fixed Bloch-sphcrc 
radius. 

In order to determine the actual squeezing limit for a 
fixed total number it is, moreover, not sufficient to take 
the general upper limit (AS2) 2 < S 2 = N 2 /4 and infer, 
from the uncertainty relation © for the variances of the 
angular momenta, the lower limit (AS3) 2 > (Si} 2 /N 2 , 
which vanishes at (Si) = and becomes as large as 1/4 
at (Si) = N/2. 

Since the variance of S2 depends nontrivially on the 
mean value of Si it is rather necessary to take into ac- 
count the constraint given by the fixed total particle num- 
ber N = 25. Using 



(S 2 ) =5(5 + 1) 
one finds, from ([5]), the inequality 

26i7 + 1 < 0, 
S(S + 1)-(S 2 ) 



S31 



with 7 = 



I (Si 



(10) 

(11) 
(12) 



This contains, besides the average, also the fluctuations 
of the "phase operator" Si . In general the above inequal- 
ity returns a lower bound £3i, m in < £31 and an upper 
bound £ 3 i,max > 61 for £31- 

We show the resulting bounds on £3 = 4(An) 2 /iV = 
4(AS 3 ) 2 /AT = 2£ 3 i|(Si)|/N as functions of the coherence 
a = 2(Si)/N in Fig. [TJ The shaded areas are allowed 
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by the inequality (fTTj) . In the medium (violet) and dark 
(blue) shaded areas the system is below the standard 
quantum limit, £3 < 1, see Eq. ([5]). 

For the special case that (Si) = S = N/2, i.e., a = 1, 
one has (S 2 ) = S 2 such that the inequality (fTTj) be- 
comes an equation fixing the variance to the unique value 
(AS3) 2 = N/A which is both a lower and an upper 
bound. Note that this bound is enhanced by a factor N 
as compared to the naive limit derived from the uncer- 
tainty relations above. For |(Si)| < N/2 one finds, since 
(ST) < 5(5+1), that 7 > 1. Moreover, (Sf) > (Si) 2 , 
and the minimum (Sf) = (Si) 2 is realized, for any (Si), 
by the eigenstates of Si. Hence, one can replace (5 2 ) by 
(Si) 2 in (fT2"j) . such that the inequality (fTT| defines the 
lower and upper bounds to £31 for any given value of the 
average (Si). For |(Si)| — > 0, one finds a lower bound, to 
quadratic approximation in (Si), of 



(AS3) 2 > 



(Si 



45(5 + 1) 



(13) 



This limit is by a factor of 5/(5 + 1) lower than the 
naive bound mentioned above, a difference which only 
disappears in the limit of large particle numbers. Also the 
maximally allowed (A53) 2 is given by 5(5 + 1) instead 
of 5 2 . 

The angular momentum representation ([3]) of the 
Hamiltonian shows that, for a given set of parameters 
U and J, each of the energy eigenstates corresponds to 
a point in the |(Si)|-(A53) 2 plane. In Fig. [2] we show 
these points for J = 1 Hz and four different values of 
U (black dots). Those values are, from top to bottom, 
U/ J = 0, 1, 10, and 100. The dotted lines serve to guide 
the eye between the values for subsequent states of the 
same Hamiltonian. This shows that states with lower en- 
ergies have larger |(Si)| and smaller (A53) 2 than states 
with higher energy. We find that the ground-state val- 
ues are in accordance with the limit set by the Heisen- 
berg uncertainty relation (jTTJ) but correspond with this 
limit only for coherences a = and a = 1. Similarly, 
the eigenstates with the highest energies, do not give the 
highest possible number fluctuations (A53) 2 . The only 
exception are the cases (Si) = ±5, i.e., U = 0, where the 
lower and upper Heisenberg limits meet at £3 = 1, and 
the case (Si) = 0. Hence, the diagram in Fig. shows 
that adiabatic changes of the parameters cannot drive the 
system, starting in some classical (diagonal) mixture of 
energy eigenstates, into the maximally squeezed state, for 
any value of the coherence a other than or ±1. We em- 
phasise that maximum squeezing allowed by Heisenberg's 
uncertainty relation is only possible in a non-equilibrium 
procedure. 

Spin-squeezed states with large total angular momen- 
tum quantum number 5 have been suggested as means 
for increasing the precision of interferometric and metrol- 
ogy measurements beyond the standard quantum limit 
[lHj]. One thereby uses, e.g., the squeezing of the un- 



certainty ellipsoid around the mean spin vector (S). If 
the ellipsoid is squeezed perpendicular to the spin direc- 
tion along a direction a this increases the measurement 
sensitivity of an angle 9 of rotation of (S) about an axis 
p perpendicular to the spin and the squeezing direction, 
(S) ■ p = a ■ p = 0. The resolution of the angle 9 is pro- 
portional to the variance (ASo-) 2 of the spin vector along 
the squeezing direction a, A9 = (AS (T ) 2 /|(S)|. This 
needs to be compared to the angular noise in the angu- 
lar momentum coherent states, 1/VN = l/\/25. Hence, 
the squeezing parameter measuring the sensitivity of the 
squeezed states considered before ((§2) = (S3) = 0) un- 
der rotations around the 2-axis reads 



N(AS 3 ) 2 

(Si) 2 + (s 2 y< 



(14) 



The area allowed by £r < 1 is indicated in Fig.[TJby dark 
shading. 

In summary, the fluctuations of the spin in one di- 
rection have to be reduced below shot noise ((AS3) 2 < 
5/2), and the spin polarization in the orthogonal plane, 
(Si) 2 + (S2) 2 , has to be large enough to maintain the 
sensitivity of the interferometer. The precision of such a 
quantum-enhanced measurement is £r/\/N [23[, whereas 
the standard quantum limit set by shot noise is 1/ y/N . 



III. DYNAMICS OF THE PRODUCTION OF 
SQUEEZED STATES 

Tuning the barrier height can be employed to produce 
many-body states with squeezing in the relative number 
difference of particles in the two wells. Such squeezing 
was observed [lH in the distribution of atoms counted 
after high-resolution imaging of the atom cloud in subse- 
quent runs of the experiment. In each such run, a con- 
densate was prepared in the potential (|Al|) with a low 
barrier height allowing the atoms to be delocalised across 
the almost flat potential floor. Then, the potential bar- 
rier was slowly raised, allowing for an almost adiabatic 
adjustment of the system's state to the modified exter- 
nal conditions. As anticipated, the finally strong barrier 
was observed to suppress the fluctuations (An) 2 in the 
particle number difference below the estimated classical 
variance (An) 2 /4 ~ N = n\ + ri2- This suppression was 
interpreted [18[ as due to the squeezing that manifests 
itself theoretically in the approximate low-energy many- 
body eigenstates. As we describe in more detail in the 
following, the near-adiabatic ramp-up of the barrier al- 
lows the atoms essentially to remain in the initially pop- 
ulated states. The change of parameters in the Hamil- 
tonian during the ramp deforms the populated states, 
causing adiabatic cooling of the system and squeezing, 
i.e., reduced relative number fluctuations which reflect 
the localisation of the particles in either of the wells. In 
a second stage when the populated levels successively be- 
come pairwise quasi-degenerate, isothermal evolution is 
observed during which the squeezing remains stationary. 



Besides a description of the dynamics observed in the 
experiment [Hj], we characterize the common properties 
as well as the differences between quantum and classi- 
cal statistical many-body evolution of the system. To 
this end we study the dynamical evolution both by di- 
rect integration of the von Neumann equation as well 
as by simulation in terms of a classical field equation of 
motion derived from the classical Hamiltonian function 
corresponding to the operator (fT|). 



A. Quantum evolution 

We consider a system whose dynamics is described by 
the Hamiltonian (fTJ). At the initial time t = to the gas is 
assumed to be evenly distributed among the wells of the 
potential, such that (S3) = at all times t > to- Its ini- 
tial state is described by a canonical density matrix with 
a given temperature To, with the spectrum determined 
by diagonalizing ([1]) for a given set of initial parame- 
ters J (to) and U(to) in the Rabi regime as discussed in 
Sect. El 

1 N 

P(to) = -^e- B '(*«'/ fc - T »|T 4 (to)>(^(^)|, (15) 

i=0 

with Z = ^2 i exp(—Ei(to)/ksTo) and Boltzmann's con- 
stant ks- The time evolution of the Bose-Hubbard pa- 
rameters J(t) and U(t) is determined by solving the 
stationary Gross-Pitaevskii equation in the double-well 
potential for a set of times during a linear ramp-up of 
the barrier Vo(t) = Vo(to) + vo(t — to) with vo/h = 
(2tt)2570Hz. Fig. [M] shows the evolution of U and J 
during the ramp up of the barrier, starting at to = 0, 
with the time given in units of the duration of the ramp 

^max • 

Decreasing in this way the tunneling parameter J and 
at the same time preserving or increasing the local inter- 
actions U drives the system from the Rabi through the 
Josephson into the Fock regime, see Sect. IA 21 We com- 
pute the corresponding time evolution of the system by 
solving the von Neumann equation for the density opera- 
tor and study the evolution of the squeezing parameters 
described in Sect. IA 31 in particular of the number vari- 
ance (An) 2 = (AS3) 2 and the coherence parameter or 
relative phase a = (S\). For the purpose of comparing 
with experimental data we fit the initial temperature To 
such that the initial number variance and coherence fit 
the experimentally determined values. 

Fig. [3] shows the evolution of a system of N = 100 
atoms in the a-(An) 2 plane as a (blue) solid line, under 
a t max = 0.16 s ramp-up of J(t)/U(t) given in Fig. [TU 
The system starts in a state with (An(t )) 2 = 37.5 and 
ot(to) = 0.994 which is obtained for an initial tempera- 
ture T = 20 nK. The (red) short-dashed line shows the 
corresponding evolution of (a(t),(An(t)) 2 ) under an adi- 
abatic change of the parameters U and J, i.e., for vo — > 0. 
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FIG. 3: (Color online) Evolution of the coherence a and the 
number variance (An) 2 for a gas of N = 100 atoms ((blue) 
solid line) under a slow ramp-up of U(t)/J(t) over a time 
tmax = 0.16 s as given in Fig. [Til The system starts in a state 
with (An(to)) 2 = 37.5 and a(to) = 0.994 which is obtained 
for an initial temperature To = 20 nK. A range of time points 
is indicated by black dots. They are spaced by 0.1 t max and 
the first completely distinct point is for t — 0.3 t max . The 
evolution under an adiabatic change of the parameters U and 
J, i.e., for vq — > is shown by the (red) dotted line. An 
isotherm for T — 0.17 nK is shown as a dashed line. See 
Fig.[T]and Sect. lA 3l for the definition of the differently shaded 
areas allowed by the constraint Heisenberg, squeezing, and 
metrology gain limits. 



The areas allowed by the constraint Heisenberg, squeez- 
ing, and metrology gain limits, respectively, are shaded 
differently as discussed in Sect. I A 31 and Fig.[T] 

Comparing the nonadiabatic with the adiabatic evolu- 
tion, as well as with the isothermal line drawn in Fig. [3] 
as a dashed line we find that the system crosses over from 
an adiabatic to an isothermal evolution. The transition 
occurs as soon as the initially occupied levels cross over 
from the linear to the quadratic part of the spectrum 
shown for different parameters U/J in Fig. [T5l 

Let us discuss this crossover in more detail. We con- 
sider the evolution of the density matrix in the energy 
eigenbasis. Fig. [4] shows the diagonal elements of the 
density matrix, (Ei\p(t)\Ej) in the basis of the energy 
eigenstates (|A7|) of the Hamiltonian (JXJ) , corresponding 
to the lowest energies Ei, i = 1, . . . , 50. Note that the 
density matrix at times t > is no longer diagonal as 
correlations have been formed. Nevertheless, the ampli- 
tudes of the off-diagonal elements are strongly reduced 
compared to the diagonal elements and the latter serve 
as a measure of the distributions of particles across the 
energy eigenstates. Fig. [4] shows these distributions for 5 
different times during the ramp-up of U(t)/ J(t) given in 
Fig. [14] on a semi-logarithmic scale. The linear distribu- 
tion at t = corresponds to a canonical density matrix 
with temperature T = 20 nK, for the linear spectrum in 
the Rabi regime, see Fig. fTSl 
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FIG. 4: Diagonal elements of the density matrix, (Ei\p(t)\Ei) 
in the basis of the energy eigenstates (|A7[) of the Hamiltonian 
([1]), for 5 different times during the ramp-up of U(t)/J(t) 
shown in Fig. 1141 Note the semi-logarithmic scale. The initial 
temperature is To = 20 nK. The inset shows the diagonal 
elements of the density matrix for t = t max as a function of 
the energy on a logarithmic scale, demonstrating the thermal 
character of the state. 



energy eigenstate during the ramp and found that the 
final state is generically far from representing a thermal 
distribution. As our results show, however, a thermal 
mixture of these initial eigenstates redistributes occupa- 
tion numbers to yield a thermal state again at the end- 
point of the ramp. The inset in Fig. [5] shows the di- 
agonal elements of the density matrix for the final time 
t = tmax in our simulations, as a function of the energy 
on a semilogarithmic scale, demonstrating the thermal 
character of the state in the levels with even index, i.e., 
the symmetric states, cf. Fig. [16J while the occupation 
of the odd-z levels remains suppressed. This suppression 
stabilizes the system against symmetry breaking to a self- 
trapped state with a non-zero mean value n. If the odd-i 
states would be equally strongly occupied, they could 
combine with the even states to yield self-trapping. 

In summary, the system changes considerably during 
the entire ramp. During the initial evolution strong adia- 
batic cooling takes place. At the crossover to the isother- 
mal evolution redistribution sets in and a "freeze out" of 
the fluctuations, fixing the system's temperature. Dur- 
ing the following evolution period, the squeezing stays 
put while the mean coherence keeps decreasing. 



B. Semiclassical evolution 



Up until around t = 0.5 t max , the density matrix 
barely changes its character indicating an adiabatic evo- 
lution. However, as long as the occupied levels are lin- 
early spaced in energy strong adiabatic cooling takes 
place, as can be inferred from the change in energy scales 
from the (red) "Rabi" to the (blue) "Josephson" spec- 
trum in Fig. [15] 

Theraftcr the system quickly leaves its adiabatic be- 
havior to develop a staircase-like shape of the energy 
distribution. Around t ~ 0.5 t max the occupied levels 
enter the quadratic regime. In this regime, the energy 
spectrum develops two-fold quasi degeneracies, starting 
in the higher levels, corresponding to symmetric and an- 
tisymmetric states, see Sect. IA 21 However, the redistri- 
bution between the states in the gradually changing spec- 
trum becomes possible because the levels at the boundary 
between the linear and quadratic parts of the spectrum 
approach each other closely. There is no redistribution 
within the quasidegenerate pairs in the quadratic part 
of the spectrum as symmetry forbids transitions between 
states even and odd in the relative particle number n. 

For times t ~ 0.8 t rnax , the energies of all but the 
largest occupied states belong to the quadratic regime. 
Note that at t = t max the occupation of the second low- 
est energy state, i = 1, is lower than the occupations of 
the ground state i = and second excited state i = 2. 
The reason for this is that the initial state p(to) is dom- 
inated by the symmetric ground state with c„.o = c_„.o 
and therefore the final state is predominantly symmetric, 
with a suppressed contribution from \E\}. 

We have studied the redistribution of an initial single- 



We will now turn to the semiclassical statistical 
description of the production of squeezed states dis- 
cussed above. This is achieved by sampling the phase- 
space probability distribution corresponding to the ini- 
tial quantum density matrix and evolving each realisa- 
tion by use of the classical equation of motion. Correla- 
tion functions at a later time are then obtained as mo- 
ments over the thus propagated probability distribution. 
The semiclassical description of the dynamics of the two- 
mode Bose-Hubbard system has recently been studied in 
Refs. [U. 

The classical dynamic equation [32j is derived from 
the classical Hamiltonian function which is obtained from 
the Hamiltonian in Eq. ([3]) by substituting the operators 
S\ = S+ + S- and S3 by the classical variables, 

S 3 ->■ n, (16) 

/ N I An 2 

S 1 ^^ 1 2 + S|cos0=-^1-— cos0, (17) 

as 



H = Un 2 - JN^jl - — cos0. (18) 

One condition for the classical description to be valid is 
that the particle number in each well is much larger than 
one, i.e., n <C N . The canonical variables are (half) the 
number difference n and the relative phase <f> between the 
two wells. From the above Hamiltonian the Josephson 



equations are obtained as 



dn 



4n 2 



(19) 



dd> An ( An 2 \ 1/2 

i = 2Un -N J { 1 -w) cos ^ (20) 

Expectation values at time t > to are obtained as mo- 
ments of the respective probability distribution 



(<D)=Af- 1 / dnd(f>P(n,<j);t)0, 



(21) 
(22) 



The probability distribution P is determined by the clas- 
sical path integral 



N= / dnd(j)P{n,(/);t). 



P{n,<t>\t)= I VnT>4>P(n,(j);t Q ) 

J to 

An 

x 8[d<p/dt-2Un--^J ( 



. ox -1/2 
iV 2 " 



An 2 

x 6[dn/dt + NJ\/l - 



sin ( 



(23) 



with the functional measures 2?n = [] T=1 dn(r), X>(^> = 
_ t d</>(r). The delta functionals evaluate the vari- 



ables n and (f> at each point in time according to the 
solution of the equations of motion, with initial values 
distributed according to P(n,,4>;to). 

For the semiclassical initial state p(to), Eq. (fl"5|). used 
in the quantum simulations in the previous section, the 
probability distribution P(n, cf>; t ) at initial time was 
determined from the Wigner function corresponding to 
p(to). See Appendix [B] for details. 

Fig. [5] shows the semiclassical evolution of a system 
of N = 100 atoms in the previously introduced a-(An) 2 
plane as a (purple) dotted line, under the ramp-up of 
J(t)/U(t) given in Fig. [Ml The initial classical distribu- 
tion was calculated from the density matrix (|15p for an 
initial temperature To = 20 nK using the Wigner func- 
tion (|B7[) . This distribution is then evolved in time using 
Eqs. (f23"l) . The (blue) solid line shows the correspond- 
ing quantum evolution as discussed in Sect. MI Al and 
Fig. [3] for the same parameters U(t) and J{t), the same 
initial temperature and the same evolution time t max . 
Obviously the semiclassical and the quantum evolutions 
are nearly identical. This indicates that dynamics of the 
production of squeezing is essentially a classical process. 
Nonetheless, the precise shape of the initial probability 
distribution as derived from the Wigner function corre- 
sponding to the state (|15|) plays a role as is shown in 
more detail in Fig. [5] for the same initial temperature as 
in Fig. [5] We plot the difference between the £3 obtained 
from the semiclassical and the quantum evolutions, nor- 
malized to the quantum result, for three different ini- 
tial phase-space distributions: An ellipsoidal Gaussian 
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FIG. 5: (Color online) Semiclassical evolution ((purple) 
dotted line) of the coherence a and the number variance 
(An) 2 = N£s/A for a gas of TV = 100 atoms under a slow 
ramp-up (t ma x = 0.16 s) of U(t)/J(i) as given in Fig.[T4] The 
(blue) solid line shows the corresponding quantum evolution, 
as discussed in Sect. IIII Al and Fig. [3] for the same evolution 
time tmax and the same time dependent parameters U and J. 
The initial temperature is 20 nK as in Fig. [3] The definition 
of the shadings was introduced in Fig. \l\ 



distribution with main-axes widths given by the num- 
ber and phase distributions derived from the state (|15|) 
(blue dashed line) , a product of the distributions of ini- 
tial relative number and coherence derived from Eq. (|15p 
(red dotted line), and a distribution as given by the full 
Wigner function for the state (fT"5|) (purple solid line). 
The frequency of the oscillations is approximately given 
by the plasma frequency in the classical potential and 
decreases with increasing U/J. We find that the solution 
derived from the Wigner function shows smaller oscilla- 
tory deviations from the exact result than that derived 
from the Gaussian distribution. The classical simula- 
tions starting from the distribution product give an even 
smaller deviation. While the distribution product gives 
the least deviations during the evolution, the initial-time 
value of £3 derived from the Wigner function is closest to 
the exact result. Its remaining deviation is due to the dis- 
crete sampling of the Wigner function. We remark that 
the oscillations take place only during the initial adia- 
batic decrease of £3, see Fig. [3l such that they remain 
mostly invisible in the comparisons shown in Fig. [5] 

The stronger oscillations of the variance of the Gaus- 
sian distribution reflect that the initial Wigner function 
which is almost entirely positive, and therefore classical- 
like, contains information about non-Gaussian correla- 
tions in the initial state. In particular, at lower temper- 
atures also quantum effects are expected to play a more 
important role which show up in the negativity of the 
Wigner function. 

In Fig. [71 we illustrate the evolution of the probability 
distribution P(n,<j>;t), Eq. ([2"3"]). for N = 100 particles, 
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FIG. 6: (Color online) Difference between the £3 obtained 
from the semiclassical and the quantum evolutions, normal- 
ized to the quantum result, for three different initial phase- 
space distributions: (Purple) solid line: Wigner function for 
the state (|15[) . (Red) dotted line: Product of distributions of 
initial relative number and coherence derived from Eq. (|15[1 . 
(Blue) dashed line: Ellipsoidal Gaussian state with widths 
given by these relative number and coherence distributions. 
See the main text for a discussion. 




FIG. 7: (Color online) Evolution of the phase-space prob- 
ability distribution P{n,(j>',t), Eq. (|23[). between t — (left 
panel) and t — 0.68 t max (right panel) as obtained by classi- 
cal simulations according to Eq. (|23[) . White color indicates 
P(n, 4>\ t) ~ while colors indicate a positive probability. The 
initial distribution is derived from a quantum gas of N = 100 
atoms at a temperature of To = 20 nK. The tilt of the fi- 
nal phase-space distributions reflects classical correlations be- 
tween the generalized position and momentum which reduce 
the squeezing in n and are due to the non-adiabatic settling 
to a thermal state, see Fig. [5] 



for the evolution shown in Fig. [5l The left picture shows 
the initial distribution corresponding to the Wigner func- 
tion (|B7|) at t = while the right panel shows P after 
t = 0.68 t max of evolution. The white areas indicate 
P(n, <t>; t) ss while colors according to the colormap in- 
dicate an increasing probability. A wide distribution in 
either direction reflects large fluctuations of the respec- 
tive observable. Given an initial distribution this would 



remain unchanged if the parameter U and J remained 
constant. Changing, however, these couplings as in the 
ramp defined by Fig.[T31 i.e., decreasing the ^-dependent 
potential term in the Hamiltonian (|18p , the distribution 
P varies as the finite distribution over "momenta" n leads 
to an expansion of the distribution in the widened cosine 
potential, see Eq. (|18[) . Tuning J to zero allows infi- 
nite expansion in the "position" direction (f> within the 
n-(f> phase space. Due to the initial finite distribution in 
n, however, classical correlations between n and (f> de- 
velop, tilting the large-time probability distribution with 
respect to the vertical axis as seen in the right panel of 
Fig. [7] and keeping the expansion finite. In summary, the 
evolution of the system in the experiment [lH can be un- 
derstood to a good approximation as classical squeezing 
of the phase-space distribution. 



C. Quantum statistical squeezing 

In the evolution of a weakly interacting quantum gas, 
quantum fluctuations to leading order enter through 
zero-point fluctuations in the initial state, i.e., they char- 
acterize the scattering into empty and out of nearly 
empty modes and play little role in the scattering in and 
out of strongly occupied modes. In the language of the 
path integral this means that the full quantum evolution 
is, to a good approximation, given by a classical Liou- 
villean propagation of the initial-time Wigner function 
which accounts for quantum fluctuations in the initial 
state, see, e.g., (3^,[3J]. Strong interactions have the po- 
tential to alter this semi-classical evolution considerably. 
However, as is illustrated by our above results, quantum 
fluctuations also in this case have only little effect if all 
available modes are strongly occupied during the evolu- 
tion. 

To see distinct effects of quantum fluctuations arising 
during the evolution requires a larger number of degrees 
of freedom part of which should stay weakly occupied. 
Alternatively, one needs to measure observablcs with a 
resolution at the few-particle level. In the two-well po- 
tential the effect of quantum fluctuations may be seen 
by introducing a tilt such that on the average only a few 
particles occupy one of the wells. 

In the following we stay with equal populations in the 
two modes but consider the variance of the relative parti- 
cle number at very low temperatures, where fluctuations 
on the order of a few atoms become relevant. At the low 
temperatures to be considered the initial state is domi- 
nated by the ground state of the Hamiltonian. Driving 
the system into the Fock regime allows to move to the 
lower left corner in the graph shown in Fig. [T] For an 
even total number of atoms the relative number fluctua- 
tions between the modes in this regime are strongly re- 
duced while the undefined phase allows interferences on 
the surface of the Bloch sphere. 

Fig. [5] compares the evolution of the exact quantum 
((blue) solid line) and semiclassical ((purple) solid line) 
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FIG. 8: (Color online) Evolution of the coherence a and 
the squeezing parameter £3 related to the number variance 
through (An) 2 = N^/4 for a gas of TV = 100 atoms un- 
der a very slow ramp-up of U(t)/J(t) as given in Eq. (|24l) . 
The quantum evolution is represented by the (blue) solid 
line while the classical statistical evolution is the (purple) 
dotted line. The system starts in the ground state of the 
Bose-Hubbard Hamiltonian for NU/J(to) — 0.16. This corre- 
sponds to (An(£ )) 2 = 24 and a(t ) = 0.993. The total evolu- 
tion time is i max = 10.18 s. The evolution under an adiabatic 
change of the parameters U and J, i.e., for 1/r — > is shown 
by the (red) short-dashed line. See Fig.[T]and Sect. lA 3l for the 
definition of the areas allowed by the constraint Heisenberg 
and metrology gain limits, distinguished by different shading. 



evolutions of a system starting in the ground state of 
the Hamiltonian with the ratio of chemical potential over 
tunneling rate being NU(t )/J{t ) ~ 0.16, for N = 100 
atoms. This corresponds to (An(to)) 2 = 24 and a(to) = 
0.993. We chose U = const, and an exponential ramp of 

J(t)/UN = J(t )/UNexp{-t/T} (24) 

with J(to)/UN = 6 and r = 0.55 s over the period of t < 
70 s. For a fully adiabatic change of the parameters U and 
J, i.e., for l/r-)0, the evolution would follow the (red) 
short-dashed line which corresponds to the dependence 
of (An) 2 on a in the ground state shown as a (red) solid 
line in Fig. 

Fig. [8] exhibits significant differences as compared to 
the semi-classical evolution shown in Fig. [5] which can be 
expressed both in terms of (An) 2 and a. With respect to 
atom number fluctuations, the classical curve approaches 
the value (An) 2 ~ 0.2 and cannot follow the quantum 
one below this limit. Following the adiabatic evolution 
for increasing initial temperatures T, we find that the 
quantum regime (An) 2 < 0.2 can only be reached for 
T < 20nK~ 10 UN. With respect to the coherence 
parameter a, the classical curve simply approaches zero 
indicating a distribution on the Bloch sphere symmet- 
ric under Si —Si. In this limit, the classical phase 



FIG. 9: (Color online) Evolution of the Wigner function 
for an inverse ramp speed r = 0.55 s, corresponding to the 
quantum evolution shown in Figs. [8land l 101 Colors encode the 
value of W(n, <j>) x 10 3 , see Eq. l(B7|. with n = (\a\ 2 - \/3\ 2 )/2 
and (j> = 2arg(a). W is negative in the blue areas. The initial 
distribution shown in the upper left panel is derived from the 
ground state of the Hamiltonian with J(to)/NU = 6 for TV = 
100 atoms. The upper right panel shows the Wigner function 
for the number-squeezed state reached at t/tmax = 0.5, with 
£3(0.5 t max ) ~ 3 ■ 10" 2 , a(0.5 tmax) ~ 0.86. In the lower 
panels we show the Wigner function at t/tmax = 0.75 (left) 
and t/tmax = 0.9 (right), i.e. for maximally negative £3 ~ 
2 ■ 10~ 5 at a ~ —0.02, and maximally negative a ~ —0.1 at 
£3 ~ 2 ■ 10 -4 , respectively (cf. Fig. [8J . The contour lines help 
to show the position of the maxima in \<j>\ = and \<j)\ = n, 
respectively. 



space distribution wraps around the equator of the Bloch 
sphere and reflects a completely undetermined relative 
phase between the wells. 

In contrast, the quantum evolution yields a coherence 
which oscillates around zero, corresponding to an asym- 
metric distribution of phases. This distribution is ex- 
hibited by the Wigner function which we show, for the 
time when a is most negative, in Fig. [9] Due to interfer- 
ences the Wigner function starts to oscillate as soon as 
the phase distribution fully wraps around the equator. 
This leads to oscillations both in (An) 2 and a during 
the near-adiabatic time evolution shown in Fig. [8l The 
oscillations in (An) 2 are damped and equilibrate at a 
value close to the point where a reached zero for the 
first time. In the lower panels of Fig. [5] the Wigner func- 
tion is shown for the evolution times t/t max = 0.75 (left) 
and t/t max = 0.9 (right), i.e. for maximally negative 
£3 ~ 2 • 10~ 5 at a ~ —0.02, and maximally negative 
a ~ -0.1 at £ 3 ~ 2 • 10" 4 , respectively (cf. Fig. |SJ). 
This shows that negative a arises from a maximum of 
the Wigner function at n = 0, \<j>\ = 7T. Moreover, in 
this pure quantum regime of the evolution, W becomes 
negative, for all <j), periodically in n, for |n| > 5. The first 
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FIG. 10: (Color online) Evolution of the squeezing parameter 
£3 related to the number variance through (An) 2 = N^/A 
for a gas of N = 100 atoms under a very slow ramp-up of 
U(t)/J(t) as given in Eq. (|24[1 . The quantum evolution is 
represented by the (blue) solid line while the classical statis- 
tical evolution is the (purple) short-dashed line. The initial 
state and all other parameters are chosen as in Fig. [8] The 
evolution under an adiabatic change of the parameters U and 
J is shown by the (red) long-dashed line. 




CO 

10T 4 r - 

10 -5 f 

0.2 0.4 0.6 0.8 1 

Normalised time t/t max 

FIG. 11: (Color online) Evolution of the squeezing parameter 
£3 related to the number variance through (An) 2 = N^/A 
for a gas of N = 101 atoms under a very slow ramp-up of 
U(t)/J(t) as given in Eq. ((24}. The quantum evolution is 
represented by the (blue) solid line while the classical statis- 
tical evolution is the (purple) short-dashed line. All other 
parameters are chosen as in Fig. \W\ The adiabatic evolution 
is again shown by the (red) long-dashed line. 



minima at \n\ ~ 6 are at the center of the grey zones. 

As was shown in Ref. [3o| the relative number vari- 
ance (An) 2 in the ground state as a function of the ratio 
NU/4J of chemical potential over tunneling coupling, 
has the approximate value (5lj 

£3 = A(An) 2 /N = (1 + NU/AJ)~ 1/2 , (25) 

in the Josephson regime where the variance is distinctly 
smaller than the classical limit, ^3 < 1, but sufficiently 
larger than 1, i.e., where iV -1 < U/J < N, sec Sect. IA 21 
Fig. [10] shows the dependence of £3 = 4(An) 2 /N on 
t, on a semi-logarithmic scale, as above for the expo- 
nential ramp (|24|) with r = 0.55 s. The (red) dashed 
line corresponds to adiabatic evolution (r — > 00) while 
the (blue) solid and (purple) short-dashed lines represent 
the near-adiabatic quantum and semiclassical evolutions 
for t = 0.55 s, respectively. The adiabatic curve shows 
three regimes in each of which log £3 varies linearly with 
time which correspond to the Rabi, Josephson, and Fock 
regimes. The full quantum evolution as shown in Fig. [10] 
exhibits an oscillatory behavior of £3, during which the 
squeezing parameter undershoots the adiabatic ground- 
state curve and eventually settles to a finite value smaller 
than that in the classical limit. Although £3 falls below 
the value it can reach in an adiabatic ramp, this does 
not contradict the Hciscnbcrg limit as is clearly seen in 
Fig.H 

As a contrast we show, in Fig. [11] that squeezing below 
the classical limit is not possible for systems with an odd 
total number N of particles. In this case, the semiclas- 
sical evolution perfectly describes the dynamics of the 



formation of squeezing correlations. To see the difference 
to the case of an even number of particles one trivially 
has to measure the particle number to better than half a 
particle. 

Taking into account the time evolution (|24p . our data 
confirms the approximate expression ([23]) in the Joseph- 
son regime. In the Fock regime the dependence of £3 on 
J is approximately given by J35| 

£3 = A{An) 2 /N = ANJ 2 /U 2 , (26) 

and the transition between the Josephson and Fock 
regimes occurs where J/U ~ 2~ 4 / 3 N~ 1 . 

Considering the non-adiabatic quantum evolution for 
J{t)/UN = (J(t )/UN)exp{-t/T} we find that &(t) 
follows the ground-state dependence as long as the ramp 
rate is smaller than the Josephson frequency (|A5|) . 1/r <C 
Lu p . As soon as the decreasing frequency uj p falls below 
1/r the squeezing parameter £3 is frozen out at the ap- 
proximate value [35j 



£3 = V 1 + (ttNUt/2) 2 - ttNUt/2. (27) 

Fig. [12] shows this limit as a function of the inverse ramp 
rate r. The solid line gives the analytical formula (|27|) 
while the (blue) crosses and (purple) circles correspond 
to the asymptotic values for £3 determined from a set of 
our quantum and semiclassical evolutions, respectively. 
Clearly, the semiclassical evolution cannot enter the Fock 
regime. 
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FIG. 12: (Color online) Final value of ^(t max ) obtained after 
an evolution time £ max = 18 At vs ramp speed r. The solid 
line corresponds to the analytical result in Eq. (|27|) valid in the 
Rabi and Josephson regimes. The (blue) crosses and (purple) 
circles are the corresponding final values obtained from a full 
quantum evolution and semiclassical evolution of the initial 
state, respectively. In the Fock regime which can be reached 
for ramp rates 1/r < 27rs~ 1 the squeezing becomes stronger 
than the classical limit £3 ~ 6 • 10 -3 . 



IV. CONCLUSION 

We have studied in detail the production of squeezed 
states in an ultracold Bose gas in a double- well trap. The 
trapping parameters are chosen as in the experiment (l8j | 
such that at the temperatures considered, the system can 
be described by a two-site Bose-Hubbard Hamiltonian. 
Following the experimental procedure, the gas is initially 
confined by a double-well trap with very weak interwell 
barrier such that free tunneling is possible between the 
sites. We studied the time evolution of the system un- 
der a slow but non-adiabatic ramp-up of the barrier, in 
particular with respect to the change in the variance of 
the particle number difference between the sites and the 
coherence which is related to the expectation value of 
the relative phase. In this way a many-body state with 
squeezing in the particle number difference, i.e., reduced 
variance of this observable at the expense of the variance 
of the relative phase, is prepared. Our results confirm 
that the squeezing attainable with a finite barrier ramp- 
up speed is limited to a value depending on the initial 
temperature and the ramp speed. This dependence is 
determined by the spectrum of the model Hamiltonian 
in which the low-energy states become quasi-degenerate 
below a certain ratio of the tunneling rate over the on-site 
energy. Once the tunneling is sufficiently suppressed such 
that the two lowest states are separated by a frequency 
on the order of the inverse ramp rate, the squeezing sat- 
urates. 

We have formulated the model and dynamic equations 
in terms of Bloch angular momentum operators and their 



correlation functions to obtain a pictorial description of 
the underlying dynamics and exhibit the connection to 
spin squeezing. Beyond a qualitative understanding of 
the experimental data of Ref. [IH our focus was set on 
the distinction between quantum and classical statisti- 
cal fluctuations. For this, classical statistical simulations 
were conducted and compared to the full quantum evolu- 
tion and description in terms of the Wigncr function. Our 
results show that within the parameter regime realized in 
the experiment, the production of squeezing is an entirely 
classical process. Squeezing below the classical limit is 
possible, however, due to the low degree of degeneracy 
in the system, only at significantly lower temperatures 
than in experiment. The detection of such squeezing and 
quantum correlations requires the measurement of parti- 
cle number at the single-particle level. Crucial differences 
arise for systems with an even total particle number as 
compared to such with an odd number. The results are 
readily applicable to other realisations which can be de- 
scribed by the model employed. We have shown that in 
the regime where quantum fluctuations become relevant, 
maximum squeezing is in general not achievable starting 
from mixtures of energy eigenstates and invoking adia- 
batic parameter changes. We emphasize that our results 
may be particularly interesting for mesoscopic dynamics 
experiments. 
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Appendix A: Bose-Einstein condensate in a 
double-well potential 

We consider a Bose condensate trapped in a double- 
well potential with two energetically degenerate minima 
separated by a barrier of variable height, allowing for 
an adjustable tunneling rate between the wells. In the 
experiment (l8j such a potential was formed optically 
by countcrpropagating laser waves creating the superpo- 
sition of standing waves with two different frequencies. 
Near the trapping minima the potential can be approxi- 
mately described as 

V eKt (x) = V cos(fcr) + \-mup-x 1 (Al) 
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x I/H FIG. 14: Time evolution of the parameters U(t) and J(t) of 

the Bose-Hubbard Hamiltonian as realized in the experiment 
FIG. 13: Trapping potential Kxt realised in the experiment This evolution drives the system from the Rabi to the 

pi| for two different barrier heights before and after an adi- Fock regime, 
abatically conducted ramp-up. 



Fig. fl3l shows the potential form for two different barrier 
heights realized in the experiment fl8l |. 

In the following we introduce the many-body Hamil- 
tonian and review its properties in the parameter regime 
chosen in the experiment [HI . We choose the description 
in terms of a two-mode model which is applicable at the 
low temperatures encountered in the experiment. For 
this model, the Hamiltonian can be diagonalised exactly 
and the number as well as the conjugate phase distribu- 
tion be studied accordingly. 



1. Hamiltonian 

As in [IH, we consider a system of identical bosons, 
of mass to, trapped in the one-dimensional potential 
|4yt (%) i Eq. (|Alj) . It is described by the Hamiltonian 

H = jdx (&(x) + Kxt(z)) *(*) 



(A2) 



Here the binary interactions between the particles are 
modeled by an effective contact potential Vi n t(x) = gS(x), 
with the coupling constant being proportional to the s- 
wave scattering length a, g = Aira/m. 

Within the single-particle sector, the ground and first 
excited states of the Hamiltonian (|A2|) . in the absence 
of a barrier (Vo = 0), are the corresponding harmonic 
oscillator cigenstatcs. Raising the barrier adiabatically 
these states are transformed into the lowest two energy 
cigenstatcs ipo(x), ipi(x) in the double well, both with 
amplitudes peaked within the two wells. As for the 
harmonic oscillator states, the ground (excited) state 
is (anti-)symmetric under the reflection r — > —r (see 
Fig. [LI. 

In the experiment [l8j], the frequency of the harmonic 
trap was chosen uj ~ (27r)5 Hz while the amplitude of the 



barrier was varied in the range between Vo ~ (27r)430 Hz 
and Vo ~ (27r)3kHz. Fig. H"3l shows the potential in these 
limits. 

For the maximum temperatures reached in [l8| , we can 
restrict ourselves to a description in terms of a two-site 
Bose-Hubbard Hamiltonian 



H 



-J(a\a2 



V 



(A3) 



where the operators di and d\ obey the standard bosonic 
commutation relations ([2]) and where J is the tunnel- 
ing and U is the onsite interaction parameter. Solv- 
ing the three-dimensional Gross-Pitaevskii equation for 
the trapping potential realized in [l8[ one finds the pa- 
rameters U and J as functions of the barrier height be- 
tween the wells. They are shown, for the potential dur- 
ing the near-adiabatic ramp, in Fig. [14] For potentials 
between the limiting cases depicted in Fig. Q2] the mini- 
mum width of the lowest band of states is obtained in the 
limit of large U/J. As will become clear in the follow- 
ing subsection, the width in this limit is approximately 
UN 2 /A ~ (27r)1.2kHz, for a total number of particles 
N = 100. To estimate the validity of the single-band 
model this needs to be compared to the typical tempera- 
tures in the experiment which are on the order of 10 2 nK, 
i.e., 10" 7 Hz. 



2. Energy spectrum 

In the following we limit ourselves to the canonical en- 
semble, with a fixed total number of particles N. A con- 
venient basis to express the energy eigenstates of (|A3|) in 
are the Fock number eigenstates of h = {a\a\ — a\a,2)/2, 

h\N/2 + n, N/2 -n) = n\N/2 + n, N/2 - ri), (A4) 

with -N/2 < n < N/2. The dimension of the Hilbert 
space is N+l. In practice, this allows us to study the en- 
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Eigenstate index i 

FIG. 15: Spectrum of the Bose-Hubbard Hamiltonian in dif- 
ferent regimes. The energies are given in units of U/s, with 
s a scale factor given in the legend. All energies are shifted 
such that Eo = 0. 
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FIG. 16: Distribution of the relative particle number n = 
(m — n2)/2 in three different eigenstates of the Hamiltonian 
in the Josephson regime, (7/ J = 1, for N = 100 particles. For 
the states with odd i, the coefficients are odd under n — > —n, 
and vice- versa for the even-i states. 



ergy spectrum by diagonalizing numerically the Hamilto- 
nian. For the following illustrations we choose a particle 
number of N = 100. 

Properties of the Hamiltonian (|3|) and the fragmenta- 
tion of a Bose-condensate in a double- well potential have 
been discussed in detail before [36j - T40| . and we review 
here only the aspects relevant to our discussion. For a 
general review see Ref. |4l| ■ See also Refs. [42|, |43[ for 
proposals for the manipulation of number and phase cor- 
relations in condensates trapped in a double well. As far 
as the form of the spectrum is concerned, three different 
regimes can be distinguished by means of the ratio U/J 
}44| . In the Rabi regime, U/J <C A^ _1 , the system con- 
sists of N nearly independent particles. This corresponds 
to the non-interacting limit. In the Josephson regime, 
N^ 1 <C U/J <C N, atom number fluctuations are small 
and coherence is high. This is called the classical regime 
as eigenstates are described by predominantly positive 
Wigncr functions with widths near the standard quan- 
tum limit. The Fock regime, N <^ U / J, is dominated by 
the interaction energy U, thus the atom number in each 
well is well defined. Reduced number fluctuations and 
other non-classical effects appear in the Fock regime. 

In the Rabi regime, the spectrum consists of a series 
of equally spaced states with a level spacing of uj p where 
the plasma frequency uj p is given by 

uj p = ^/2J{NU + 2J) (A5) 

Increasing U / J beyond 1 /N introduces an approximately 
quadratic part to the spectrum for energies E > 2NJ, 
for which 

E. t ~NJ+ji 2 . (A6) 

In the Fock regime only quasi-degenerate pairs of states 
remain, forming an approximately quadratic spectrum. 



It was shown in Ref. [45( that the splitting between these 
quasi-degenerate states vanishes with 1/Nl. 

Fig. ll6l shows the occupation number distribution |c n | 2 
of three different energy eigenstates \E/) in the Fock basis 
|ra) = \N/2 + n, N/2 - n) defined in Eq. flA4}, 

N/2 

\Ei) = ( A? ) 

n=-JV/2 

We find that the coefficients of the ground state are well 
approximated by a Gaussian centered around n = 0, re- 
flecting its semiclassical nature. 

The first exited state has a wider distribution around 
the same mean value. Since it must be an antisymmetric 
state, one has ci.o = 0. The higher excited states show 
an increasingly wider distribution around i = with an 
increasing number of "nodes" in c n ,i- 

In the Fock regime, U/J > N, all states are quasi 
twofold degenerate. In this regime, there is no longer 
any contribution near n = 0. A gap opens up in the 
distribution |c n .;| 2 which increases towards larger i. In 
the limit U/J — > oo the highest excited state, 

\E N+1 ) = ^=(\N,0) + \0,N)). (A8) 

is cat like. Note, however, that excitation of the double- 
well system up to \En+i), e.g., at higher temperatures, 
easily exceeds the validity of the single-band approxima- 
tion. 



3. Angular- momentum representation 

In many-body quantum physics, number and phase 
are a particular type of conjugate variables which ful- 
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fill a Hciscnbcrg-type uncertainty relation. We discuss, 
in this subsection, the representation of the Hamiltonian 
in terms of number and phase operators and review dif- 
ferent measures of squeezing. The uncertainty relation 
allows us to derive fundamental limits for the respective 
squeezing parameters. 

The relation between number and phase can be most 
conveniently visualised on the Bloch sphere, i.e. in terms 
of the Schwingcr representation [46[ of the SU(2) sym- 
metry group, 



Si- 



1 



2 



k = 1,2,3, 



(A9) 



where cr k is the Pauli /c-matrix. As a consequence of 
the commutation relations (|2|) , the operators Si form the 
fundamental representation of the angular momentum al- 
gebra, with 



[Sk, Si] = ieklmS-n 



(A10) 



£fcZm being the total antisymmetric tensor of rank 3. 
These commutation relations give rise to a set of un- 
certainties relations that determine lower bounds for the 
fluctuations of the observables, 



(AS k ) 2 (AS l 



> -^\£klm{Sm)Y 



W) 



with (A5fc) 2 = ((Sk - (Sk)) 2 )- In this representation 
the two-mode states with fixed total particle number N 
can be written as angular momentum states \S, S3) with 
S = N/2, S 3 = n. The Hamiltonian (fAl]) in this repre- 
sentation is given in Eq. ((3j) . 



The antinormally ordered characteristic function is then 
obtained by Fourier transforming the Q-function, 



C A (X,p) = j d 2 Qd^Q(a,/?)e At,, - A *V r -"*' ? . (B4) 

The symmetrically ordered characteristic function is fi- 
nally calculated using the Baker-Campbell-Hausdorff re- 
lation, leading to 



(B5) 



Cs(\^) = C A (\^)e^ 2+ ^/\ 



and from this the Wigner function by use of Eq. (|B1[) . 

Applying this procedure and using the representation 
of the density matrix in the Fock basis (|A4|) . |n) = \N/2+ 
n,N/2-n), 



\n\p\m) 



(B6) 



we arrive at 



W(a,(3) =- 



N 

£ 



Pnr 



yJn\m\(N - n)\(N - mj\ 



x Q nm (a)n N _ n w - m (/3)e- 2(|a|2+l/32|) , 



(B7) 



with 



Appendix B: Wigner function of two- mode states 

The Wigner function of a two-mode system can be cal- 
culated from its definition 



W(a,/3) 



1 



d 2 Ad 2 uC s (A,^)e 



a\* —a* \ jSfi* —{3* fi 



(Bl) 

where Cs is the symmetrically ordered characteristic 
function, which one obtains as follows [471 ] . See also 
Ref. [HI for a detailed discussion of two-mode Wigner 
functions on the Bloch sphere. First, the Q-function 



Q(a,P) = \(a,P\l>\t*,P) 



(B2) 



fc=0 ^ ' 1=0 ^ 



(B8) 



x Hk+i(2lu\a)H n+m _ {k +i){2Rca) 



where H n (x) is the nth Hcrmite polynomial. Eqs. (|B7|) 
and (|B8[) show that W is real if p is hermitian. 

As we only consider closed systems, the dependence 
of W on the absolute phase of the two modes is irrele- 
vant. We therefore need to evaluate W only for different 
relative phases 0, choosing, e.g., 



a = \a e 



| a |e 1 */ 2 , p = \fi\e~^' 2 . 



(B9) 



is determined as the expectation value of the density ma- 
trix p with respect to the two-mode coherent state 



|j8 ) =e -(MW)/2 ^ 



a 1 ^ 



i,j=0 



\i)i\j)r- (B3) 



Taking furthermore into account that the total number 
of particles N = \a\ 2 + |/3| 2 is fixed reduces the number 
of free arguments of W to two, the quantities n = (\a\ 2 — 
|/3| 2 )/2 and 4> = 2arg(a). 
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